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1.  Introduction 

Let  % = {Xi, . . . , } be  a point  set  of  independent 

identically  distributed  (/-dimensional  random  vectors 
sampled  from  the  probability  measure  and  AT  be  a 
punctured  at  the  origin  cone  in  R'',  (/  > 1 . We  define  the 

kth  layer  as 

= {X,  : # {Kx,  n = it-l  } it  = 1,2, , 

where  = x+K  is  the  translated  cone  with  vertex  in 
X E R‘^.  Intuitively,  the  /:th  layer  is  the  set  of  the  A:th 
extremes  of  in  the  direction  X.  The  prime  examples 
we  have  in  mind  are  (1)  the  Pareto-optimal  points  corre- 
sponding to  the  first  layer  in  the  direction  of  the  positive 
orthant,  and  (2)  the  total  maximum , which  may  be  con- 
sidered as  the  first  layer  in  the  direction  of  the  cone, 
complement  to  the  negative  orthant.  We  are  interested 
here  in  the  distributions  of  random  variables 

counting  the  number  of  points  in  the  /:th  layer.  These 
distributions  depend  essentially  on  both  K and  p. 

From  a more  general  viewpoint,  the  first  layer  can  be 
regarded  as  the  set  of  maximal  elements  [4]  with  respect 


to  the  binary  relation  '3i  in  defined  as  x^y  <=>x- 
y E K.  A Iternatively,  any  scale  and  translation  invariant 
binary  relation  generates  a cone  by  setting  K = {x  E R*^: 
x0lt)}  and  the  maximal  elements  are  conical  extremes. 

Two  above  cases  of  the  counting  problem  have  been 
considered  in  the  literature  under  the  assumption  that  p 
is  either  a product  of  one-dimensional  marginal  mea- 
sures or  a multivariate  normal  distribution  [2,10,1 1,12]. 
It  is  well  known,  for  example,  that  if  yii  is  a product 
measure  in  R''  then  the  average  number  of  Pareto  points 
is  of  the  order  of  (log  n while  the  probability  that  the 
multiple  maximum  exists  is  n 

In  this  paper  we  focus  on  a class  of  distributions  p 
already  studied  in  connection  with  multivariate  ex- 
treme-value theory  [8]  and  statistics  of  convex  hulls 
[1,5 ,6,9].  These  distributions  are  characterized  by  regu- 
lar variation  of  the  tail  of  the  radial  component  and 
asymptotical  independence  of  radial  and  angular  com- 
ponents. We  show  that  typically  the  V„‘*’’s  converge  in 
distribution  and  the  expectations  have  finite  limits  as 
n ^ 00.  In  the  special  case  of  slow  variation  we  calcu- 
late explicitly  the  limiting  distributions. 
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2.  Preliminaries 


and  using  the  iid  assumption  write  for  the  expectations 


We  define  a cone  as  a punctured  at  the  origin,  scale- 
invariant  Borel  set  in  R‘^,  i.e.,  0 ^ K,  tK  = K V/>0. 
Each  cone  is  uniquely  determined  by  its  spherical  base 
S+  = ^ n S , where  5 denotes  the  unit  sphere.  We  asso- 
ciate with  K also  the  spherical  set  5^  obtained  by  reflec- 
tion about  the  origin,  So  = 5\  (5+  U S_)  and  S±  = S+  H S- 
The  cone  with  spherical  base  CCS  will  be  denoted 
cone  (C). 

Set  B,=  {xER'':l|.r||<r},  = R^S,  and 

Ar,c  = cone  (C)  H 

We  fix  in  what  follows  a cone  K and  a multidimen- 
sional probability  distribution  fx  satisfying  the  following 
conditions: 

(i)  There  exists  a > 0 and  a probability  measure  p on  S 
such  that 


lim 

/— >00 


r > 0, 


(1) 


(ii)  For  all  p-continuous  CCS 


hm  =p(C), 

,->00  p,(S,) 


(2) 


(iii)  p (int  S+)  > 0,  and 


(iv)  p has  no  atom  at  the  origin. 

Consider  an  iid  sample  from  p,  ^„  = {Xi, . . . , X„  }, 
represented  in  the  polar  form  as  the  product  of  radial 
and  spherical  components:  X,  = /?,  Z,,  where  /?,  = ||X,||, 
Z,  = X//|lX,||.  The  above  conditions  on  p have  a natural 
probabilistic  interpretation.  Condition  (i)  means  that  the 
distribution  function  of  the  radial  component, 

Fir)  p (B,), 


has  a regularly  varying  tail.  Condition  (ii)  is  translated 
as 

lim  P{Z|  C ■\Ri>  r}  = p (■) 

r-^cc 


n B)  = nP{x,  e n B } = 

n ) P{X,  GB  ;X2,. . . ,X,GXx  ,;X*>i,.. . ,X„  }= 

The  following  lemmas  will  be  used  to  estimate  these 
integrals. 

Lemma  1.  There  exists  t>  0 such  that  p(X^)  > 
T(l-p(BH))/or  all  xCR\ 

Proof.  Consider  first  the  case  where  there  exists  a 
linear  isomorphism  which  maps  K onto  the  positive 
orthant.  Let  y be  the  inverse  image  of  the  vector 
(!,...,!)  under  this  isomorphism.  By  convexity, 
KyCKJor  alljrGB,. 

Condition  (iii)  allows  one  to  select  a compact  p-con- 
tinuous set  CC  S+  with  p (C)>0  . It  is  easy  to  see  that 
y G int  X,  the  sets  K^y,  s>0,  are  increasing  as  5 f 0 and 
Uj>o  Ksy  = int  X.  It  follows  that  C C Ksy  for  sufficiently 
small  s.  Furthermore,  for  small  s we  have  also 
Ai  c C Xjj,.  Indeed,  the  points  of  i4,.c  are  representable  as 
tx  , with  t > 1,  jc  G C,  thus,  by  convexity,  x G Ksy  im- 
plies tx  G K,sy  C Ksy.  Homogeneity  implies  A„s,c  C X^. 
It  follows  now  from  Eqs.  (1)  and  (2)  that 

P'iKry)  P'iArls.c)  P'iArls.c)  piPrls) 
p(B?)  ""  p(Bf)  " fji(B^,s)  p(Bf) 

-4  5“p(C),  . (4) 

From  Kj,  D X^i^  we  derive  for  sufficiently  large  ro  and 
ll^l  > ro  that 

p(X,)  ^ p(X|Hb)  1 
p(By-  p(By  2 

For  X G Bro  we  have  Kj^  D Kr^y , therefore  Eq.  (4)  along 
with  the  inclusion  K,r„y  D Kr„y,  t > 1,  implies 


and  is  to  be  interpreted  as  the  asymptotic  independence 
of  radial  and  spherical  components,  where  the  limiting 
distribution  p does  not  disappear  in  the  interior  of  5+ 
(condition  (iii)).  The  last  condition  is  not  essential  and 
assumed  for  technical  reasons. 

Given  a Borel  set  B C R*^ , we  represent  the  number 
of  the  /:th  layer  points  in  B as  the  sum  of  random  indica- 
tors 

n 

n B ) = 2 1 <x,  e n b}5 

/-I 


P'jKx) 


pi^roy  ) 


> p,iKrgy)  > 0 . 


The  assertion  follows  in  this  case  by  setting 
T=min(p(X,„,),^5“p(C)). 

For  arbitrary  X one  can  find  a smaller  cone  X'  C X, 
which  is  linearly  isomorphic  to  the  positive  orthant  and 
still  has  the  interior  of  its  spherical  base  of  positive 
p-measure.  This  is  possible  since  the  spherical  d-sim- 
plexes  build  a measure-generating  class  on  5.  It  remains 
to  note  that  p(X^)  > p.(Xj)  for  any  translation,  whence 
the  estimate  holds  in  general.  □ 
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Lemma  2.  //E  has  a limiting  value  v G [O,*) 
then  all  EK®,  k = 2,3,  • • • converge  to  this  limit  as 

Proof.  Let  m{t),t  E [0,i  ],  be  the  distribution  ftinction 
of  the  image  measure  induced  by  the  mapping 
Changing  variables  transform  Eq.  (3)  to  the 
one-dimensional  Lebesgue-Stiltjes  integral 


Similarly, 

lim  sup  (l-T(l-F(r)r' 

«-^°o  Jo 

dF(r)= 

r{y„)  r'' 

(i-T(i-t)r'df(t)-n  (i-T(i-f)r'  di 

) -^0 


EV®  ^ 


trMm(f). 


where  the  equivalence  can  be  justified  by  partial  integra- 
tion. □ 


A slight  modification  of  the  standard  Tauberian  theorem 
as  found  in  [14]  assures  that  the  limiting  value  of  this 
integral  for  i = 1 exists  iff  m(t)  is  left-differentiable  at 
r = 1 , in  which  case  the  limit  and  the  derivative  have  the 
same  value.  Applying  this  theorem  in  the  reverse  direc- 
tion one  can  easily  see  that  all  the  EV„®  ’s  must  have  the 
same  limit.  □ 

Lemma  3.  Assume  jn  ij  an  increasing  sequence  such 
that  lim  n{l-F , j > 0 , then 

lim  sup  EV®  ^ T'*,  (5) 

lim  sup  E # ) n ^ e-^T“ (6) 

n^oo 

with  T determined  by  Lemma  1 . 

Proof.  Set 


/(r)=  dF(r).  rG[0,l]. 

Jpir)  S f 

Regular  variation  of  F at  infinity  implies  readily  that  for 
all  sufficiently  distant  discontinuity  points  the  ratio 
(iump-size)/(distribution  tail)  is  close  to  zero.  It  follows 
that  (l-/(f))/(l-f)-^l  as  r T 1 (for  continuous  F this  is 
obvious  since  fit)  ==  t ).  Lemma  1,  a change  of  variables 
and  the  Tauberian  theorem  yield 

EV„®  = Rf  {l-fjLiK,)r-'dfjL(x)  ^ 

Jr*' 

n[  (l-T(l-/x(FH))"“‘djLt(x)  = 

jRrf 

[(l_7-(i-F(r))''-'dF(r)=  [ (I-T(l-t))"-‘d/(r)^T-\ 

Jr  Jo 


3.  Pareto-Tailsj  a > 0 

In  this  section  we  study  the  limiting  behaviour  of  K® 
under  the  assumption  that  the  regular  variation  index  a 
in  Eq.  ( 1 ) is  positive.  Our  plan  is  to  translate  Eqs.  ( 1 ) and 
(2)  into  the  convergence,  of  a suitably  normalized  sam- 
ple, to  a Poisson  process[6,15]  and  then  apply  a continu- 
ity argument  to  prove  also  the  convergence  of  the 
V„®’s. 

Compactify  R*^  by  adjoining  the  infinite  point  oo  and 
then  puncture  in  the  origin.  The  resulting  topological 
space,  say  is  isomorphic  to  R''  and  canonically  em- 
bedded into  its  compactification,  bounded  from  the 
origin  Borel  sets  B CW  being  relatively  compact.  We 
endow  the  space  M(R‘^)  of  Radon  measures  with  the 
vague  topology:  m„  m iff  m„  (B  )—^m  (B ) for  all  rela- 
tive compacts. 

There  exists  a sequence  of  positive  constants 
such  that  the  measures  v„{-)  £ njx{a„-)  converge 
vaguely  to  the  measure  v determined  by 

KA,c)  = r-XC),  K{°°»  = 0.  (7) 

The  limiting  measure  is  in  MiW^)  , being  infinite  on 
balls  centered  at  the  origin  as  well  as  on  the  sets 
cone  (C)  with  p(C)>0.  In  particular,  condition  (iii) 
implies  v (int  K)  = <x>.  Clearly,  visa  product  measure  in 
polar  coordinates  and  has  no  atoms. 

Let  ^ be  a Poisson  point  process  in  R"^  with  intensity 
measure  v,  and  be  the  random  element  of  MiW) 
associated  with  the  scaled  sample  af^„.  Obviously,  the 
operation  of  taking  a layer  commutes  with  rescaling: 

a>0,  therefore  the  number  of 
points  in  each  layer  remains  invariant  under  scale  trans- 
formations. One  can  expect  in  this  situation  that  V„® 
converges  in  some  sense  to  an  analogous  functional  of 
the  Poisson  process. 

Define  the  ^th  layer  of  the  Poisson  sample  as 

= {x  G R"  : f ({x})  = Li  (KJ  = k-l}, 
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Using  the  polar  representation,  x = rz,  and  homo- 
geneity we  can  write  = r~°‘4>{z),  where  ^{z)  = 
v(KJ  is  a function  of  the  spherical  argument.  Using 
Palm  probabilities  and  integrating  along  radial  rays  we 
represent  the  expectations  as 

(ft—  1)1 


Js 


I 


dp(z) 
4>iz)  ■ 


(k-1)! 


dp(z)d(-r  “)  = 


(8) 


The  resulting  integral  does  not  depend  on  as  it  is 
suggested  by  Lemma  2.  The  integration  area  can  be 
reduced  to  5\int  5_  since  cf)  is  infinite  on  int 

The  following  lemma  is  found  in  [5]. 

Lemma  4.  Let  E be  a locally  compact,  Hausdorff  and 

separable  space;  ho,  h be  a uniformly  bounded 

sequence  of  real  measurable  functions  commonly 
supported  by  a relatively  compact  set;  and  mo,  m\, . . . 
be  a sequence  of  Radon  measures  on  E such  that 
m„  A mo.  The  set  D = {x  E : 3{jr„  },  x„^x,  h„(x„) 
-h  h{x))  is  measurable  and  if  mo(D)  - 0 then 
f h„dm„^f  hodmo. 

Now  we  are  ready  to  prove  a convergence  result. 


Theorem  1.  Assume  (i)-(iv),  a > 0,  and 

vi-dK)  = 0,  (9) 

vXv{{x,y)  G R^XR'^:  (x-yj  G dK}=  0.  (10) 

Then  for  all  k = 1,2,... 

(v?,...,0-A(V">,...0,  (11) 

/2^^.  (12) 


Proof  By  Skorohod’s  theorem  we  can  find  random 
point  measures  ^„,^GM(R'')  satisfying  = 

I ^ ^and  In  ^ a.s.  Thus  to  prove  the  convergence  in 
distribution  (10)  it  suffices  to  consider  the  case  ^ 
a.s.  In  what  follows  we  fix  a typical  realization  of  ^ and 
assume  n sufficiently  large. 

Since  r-lint  K)  = ^,  ^ lays  in  the  cone  interior 

infinitely  many  points.  Select  k of  them,  say  Xi x*. 

Pick  r sufficiently  small  to  satisfy  Br  C Dj. , as 
well  as  = 0 and  also  Br  D {xi, . . . , x^ } = 0.  The 

complement  is  relatively  compact  hence  the 
processes  and  ^ have  there  a finite  number  of  points. 


say  y„,u  ...  , y„.p  and  yi, . . . , y^,  respectively.  These 
points  may  be  labeled  so  thatyn,,  — > y, , as  it  follows  from 
the  vague  convergence.  By  the  construction,  any  trans- 
lated cone  with  x E Br  contains  at  least  k points  of 
thus  Br^  5e^\^)  = (d  and  also  5,  D = 0 for 

j=l,...,k. 

For  y,  G -int  /( , the  cone  Ky.  contains  a vicinity  of  the 
origin,  where  ^ has  infinitely  many  points.  Therefore 
y,  ^ nj. , yn..  ^ n‘. , 

The  condition  shown  in  Eq.  (9)  assures  that  no  one  of 
yi, . . . ,yp  lies  on  -SK,  almost  surely. 

For  y,  ^ -cl  K,  the  shifted  cone  Ky^  is  bounded  from 
the  origin.  Therefore  there  exists  an  open  vicinity  of 
cl(U/’=  1 KyJ  which  is  still  relatively  compact,  and  hence 
contains  at  most  a finite  number  of  points  in  addition  to 
yi, . . . , yp  . By  Eq.  (10),  ^{dKy.)  = 0 a.s.  Again  the 
poinwise  convergence  implies  ^ = 

^{KyX  whence  (V„“’, . . . , V/*')  = (V‘", . . . , V<*>)  and 
thus  (11). 

Now  turn  to  the  convergence  in  mean.  It  is  enough  to 
prove  Eq.  (12)  for  the  first  layer,  A:  = 1.  It  is  easy  to  see 
that 

EV„“’  ~ j exp(-u„(KJ)du„(x)  - j ( . . . Idr-n-i- 
I (...)di^„  r>0 

Take  a point  x ^ -cl  A"  and  a sequence  x„  — > x and 
consider  the  indicator  functions  of  the  sets  K^c  and  Kx„  as 
the  h's  in  Lemma  4.  The  divergence  set  D is  dK^, 
whence  by  (9)  and  the  lemma  u„  (Kx„)  v{Kx). 

For  X G - int  K,  x„— >x  we  have  Vni^xJ  to  viKJ  - °o 
since  contains  some  fixed  vicinity  of  the  origin,for 
all  sufficiently  large  n.  Therefore  in  this  case  also 
K {K,„)-^v{Kd  = 00. 

To  make  further  use  of  Lemma  4,  consider  this 
time  the  functions  ho{x)  = expl-r'i/f.t)),  /t„(x)  = 
exp(-r’(A'^„)).  For  the  discontinuity  set  we  have  D C 
-dK  U {x  : jfdK^)  >0}.  The  assumptions  in  Eqs.  (9) 
and  (10)  imply  ifD)  = 0 , hence  for  any  r 


■’b; 


exp(-i^„(A:U)dr'„(x) e\p{-viK,))dv(x). 


Now  apply  Lemma  3 to  derive  the  estimate 


lim  sup  exp(-j^„(A:j)di^n(Jf)  ^ T''exp(-i^(fi,^)T). 

The  right-hand  side  here  tends  to  zero  as  r— >0. 

Putting  this  all  together  and  comparing  with  Eq.  (8) 
we  conclude 

lim  sup  EV‘,1'  ^ EV^'K 
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The  reverse  inequality. 


lim  inf  > EV“>, 

/I— >0C 

follows  from  the  convergence  in  distribution  □ 

Remark.  The  continuity  conditions  of  Eqs.  (9)  and 
(10)  are  actually  some  properties  of  the  spherical  mea- 
sure p.  The  first  one  trivially  translates  as  p (-BK)  = 0, 
but  we  have  not  been  able  to  find  a re-formulation  for 
the  second.  Sufficient  conditions  for  Eq.  (10)  are;  p is 
non-singular,  and  dS+  lies  in  a (n-2)-dimensional  set; 
or  is  convejt,  SK  has  no  two-dimensional  facets  and 
p(^5.)  = 0. 

Example.  Here  is  a remarkable  case  where  the  expec- 
tations are  explicitly  computed.  Consider  the  two- 
dimensional  Cauchy  distribution  specified  by  the 
density  dp,(x)  = (27r)"'(l-i-||jc||)“^'^djr,  x G R^.  The 
radial  tail  is  regularly  varying  with  a = 1 and  the  circu- 
lar measure  is  uniform,  i.e.,  dr'  (rz)  = {2tt  )“‘r“Mrdz, 
r > 0,  z £ [0,2 tt). 

Assume  first  that  K is  the  positive  quadrant.  The  ^th 
layer  are  those  X,  ’s  which  are  exceeded  by  exactly  k~  1 
points  of  Mn  in  both  components.  Integrating  yields 


ct>{z)=viK,)  = 


cos  z +sin  z -1 
2tt  sin  zcos  z 


z G {-ttI2,  tt) 


and  </)(z)  = <=°  otherwise.  Computing  the  integral  in 
Eq.  (8)  we  obtain 


lim  EV„“>=  1 +~. 

4 


The  first  layer  is  either  empty  or  just  one  point,  maxi- 
mizing both  components,  thus  this  mean  value  coincides 
with  the  limiting  probability  of  the  total  maximum. 

The  limiting  distribution  and  higher  moments  of  the 
’s  can  be,  in  principle,  expressed  in  terms  of  some 
integrals  similar  to  Eq.  (8).  These  expressions  do  not 
seem  tractable  by  analytical  methods  because  of  the 
complicated  integration  domains. 


4.  Slowly  Varying  Tails:  a = 0 

The  case  of  slowly  varying  radial  tail,  with  a = 0 in 
Eq.  (1),  is  of  special  interest.  The  above  Poisson  approx- 
imation method  does  not  work,  since  the  sample  cannot 
be  rescaled  to  provide  a non-degenerate  limit.  To  get 
around,  we  extend  here  a method  already  exploited  in 
[1],  where  the  number  of  convex  hull  extremes  of  a 
sample  under  slighly  stronger  assumptions  on  the 
distribution  has  been  studied. 

We  assume  for  technical  reasons  that  F is  continuous 
though,  in  fact,  slow  variation  is  all  that  is  needed. 

Let  . . . , X';?'  be  the  elements  of  arranged  in 
the  norm-decreasing  order,  i.e.,  ||X';1'  1|  > . . . > ||X';?^  1|. 
Set  = Rj  and  Z''*  = Z, , iff  X*''  = X- i,  j = , n. 

One  can  recognize  in  the  R ’s  the  radial  order  statistics. 
The  associated  spherical  variables,  Z'^'  ’s  will  be  called 
concomitants . Note  that  the  continuity  hypothesis  make 
the  definitions  correct  since  the  radial  components  are 
different  with  probability  one. 

Mailer  and  Resnick  [13]  proved  that  slow  variation  is 
equivalent  to 


Rl;^ 


/=  1,2, . . . . 


(13) 


For  /:  = 1 we  have  the  limiting  mean  of  the  number  of 
Pareto  points. 

Now  suppose  K is  the  complement  to  the  negative 
quadrant.  The  /:th  layer  consists  of  those  X,’s  which 
exceed  all  except  some  k- 1 sample  points  in  both  com- 
ponents. We  get 


Our  convergence  results  effectively  exploit  this  fact 
combined  with  the  asymptotic  independence  of  the  con- 
comitants shown  next. 

Let  Z''*,  Z'^’, ...  be  iid  S-valued  random  variables 
with  distribution  p. 

Lemma  5.  Assume  that  F{r)  = pi{Br)  is  continuous 
and  Eq.  (2)  holds.  Then 


(piz)  = 


cos  z -H  sin  z -Hi 
2tt  sin  z cos  z 


Z G (0,77/2) 


...,  Z^)  ^ (Z"', . . . ,Z'*')  k=\x... 


and  (/)(z)  = =«  otherwise.  Computing  the  integral  Eq.  (8) 
in  this  case,  we  obtain 


lim  EV„‘'^>=  1 
n—^oc  4 


Proof.  For  p-continuous  CCS  write  (2)  as 
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We  have 


Introduce  the  random  variables 


P{Z'^'  E Cu  = X,, . . . , Z'i'  E Q,  = 

X,}  = ^^^P{Z'i'EC„ /?'■'  = /?, 

Z^  E Q , = /?*}  = P{Z'^'  EC,,..., 

Z'*‘  E C*  ,/?!>...>  ;/?*>/?,  for  / = n } = 

j (FCrOrMyu-iriZ,)  . . . d/x(nzi)  = 

Ti  > . . . > rt 

z.  GC, z*eG 


r„'‘»  = 0;  r„<'>  = min{/  : uLi  Cl  {Xj'^", 

...,XW>  = 0},  k=  !,...,«, 

which  count  the  X,  ’s  in  the  norm-decreasing  order  until 
the  first  k layers  having  been  filled.  Clearly, 
yw  ^ ^ Denote  5“  the  product  of  infinitely 

many  spheres,  and  set 

f “»(z)  = 0;  = min{/  : (5.  Cl  {z„  ...,Zj})  = k}, 

= #{/•  : f**-"  (z)  <;•  ^ f'"(z),  z;  E S\S^}, 

where  z = (zi,  Z2, . . . ) ^ >5“  and  inf  0 = 00.  For  i S: y 
the  set 

{z  E 5”  : t/'“(z)  = y , f^iz)  = / } 


n! 


(«-/:)! 


I 


(F(r*)r‘dFc,(r,)  . . . dFc,(r,)  = 


n-\ 

k-\ 


I 


"(ri)  I 


(F(/-*)r'dFc,(A-,) . . . dFc,(n)  = 


r > Ti 


(F(n))''-MFc,(n) 


Li 

- "(ri)  I 


dFc,(r,) 

(F(n))"-‘ 


. . . dFQ.,(n.,) 

(1-Fc,  (n)) . . . (1-Fc,_,  (/-,))  dFc,(/-.)  = 


is  a finite-dimensional  cylinder  in  5”. 

Denote  Z = (Z'‘*,  Z'^’, . . . ) the  sequence  with  iid 
components  distributed  according  to  p,  T*’  = 
y(i‘)  ^ v(*>(Z).  It  follows  from  the  definition  and  condi- 
tion (iii)  of  Sec.  2 that  = 1 ,2, ...  is  a strictly 

increasing  sequence  of  finite  stopping  times  with 
respect  to  Z''',  Z'^’, . . . 

Theorem  2.  Assume  (i)-(iv),  a = 0,  p(<9S+  U <?S_) 

= 0.  for  any  A = 1 ,2, . . . 

T,!’, . . . r'*>)  ^ (v“>,  r“>, . . . , v<‘>,  py 

Proof.  Fix  integers  v,, . . . , v^;  r,, . . . , = r satisfying 

0 </,<...<  r*  and  0 < v,  for  / = 1 , . . . , A.  We 

need  to  prove  that 


"to 

(for  large  r uniformly  in  n ) 

"(rO  J (l-F(r,))*-'(F(r,)r* 

dPctil'k)  + € ^ P(C|)  . . . p(Ck)  + €, 

as  , where  we  have  used  Eq.  (14)  and  applied  an 
argument  similar  to  that  in  Lemma  3.  Asymptotically, 
the  probability  is  factorized,  whence  the  statement  □ 
To  prove  the  convergence  we  combine  in  what  follows 
the  above  lemma  and  Eq.  (13).  The  idea  is  that  the  points 
with  top  layer  ranks  have  also  small  ranks  in  the  radial 
components.  On  the  other  hand,  conical  extremality  of 
the  points  with  small  layer  ranks  is  determined  by  their, 
almost  independent,  spherical  components. 


lim  P{V<'>  = V,,  7’^;’  = ?,;  /^A}  = 

/!— >00 

P{F«  = v,,  r<'>  = r,;i^A}. 

We  endow  S'  with  the  product  measure  p'  and  the 
Euclidean  metric.  Define 

D = dS^U  SS.,  3)  = {(zi, . . . , z,)  G 5' : 

{z„...,z,}nD  =0}, 

,p^={(z„...,z,)e5':  V'«(z„...,z,)  = 

Vi,P‘\zi,...,z,)  = tr,i^k}. 

The  definition  of  M is  correct  due  to  the  cylindrical 
property.  It  is  easy  to  see  that  3 is  compact,  si\3  is 
open  and,  by  the  assumption,  p'(3)  = 0.  It  follows. 
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p'{M)  = p’{si\3}).  For  any  8 there  exists  s4s  with  the 
properties: 

sisC.M\3),  p(si}-p(siis)<8,  dist(,5^s,®)>0, 
pXd^s)=0.  (16) 

To  prove  this,  take  %={aE,5A®  :dist  {a,S(6^ 
2)}  >6},  then  is  an  open  set,  increasing  to 
,52^^  as  . We  have  p'{{,s^)\%)  < 8 for 

sufficiently  small  0.  On  the  other  hand,  C 
{a  G 5' : dist  (a,^(42f\®))  = 0),  these  sets  being  dis- 
joint for  different  9.  Hence  the  set  of  the  values  of  9 
with  p'{d‘^e)  > 0 is  at  most  countable.  Select  an  appro- 
priate 6 and  set  '^g . 

We  derive  from  Eq.  (16)  with  the  help  of  some  topo- 
logical considerations  that  for  sufficiently  small  e 

LJ  (z/ + fid  n cone  (D ) = 0 (zu  ■ ■ ■ , z,)  E:  sis- 0^) 

i ^ t 

Assume  now  that  the  compound  event 

(Z'i), . . . , Z'^)  G Sis-,  (18) 

occurs.  We  show  next  that  in  this  case 

Tp  = t,-, y«  = v,.  i=\,...,k.  (19) 

Let  Q be  an  element  of  the  finite  algebra  of  spherical 
sets  generated  by  S+  and  S_.  The  following  equivalence 
holds: 

e cone(e )«Z'i’ G g for  1 ^ < r ,i< n . (20) 

Indeed,  note  first  that  dQ  C D.  By  (17),  Z'^l}  G Q implies 
Z^,;’ -t- fie  C cone  (0).  From  Eq.  (18)  we  have  also 
X'iNfi;j^iCcone(e).  But  -X*'*  G fi^^i  thus  X''=- 
Gcone(Q).  Use  Q""  instead  of  Q to  prove  the 
reverse  implication. 

The  definition  of  si  and  Eq.  (18)  yield  Z'2’  G S+, 
Setting  g = S+  in  Eq.  (20)  we  have 
X^^  E.  K , ti<j  ^ n.  Therefore, 

{Xl''"“, . . . , X'”’}  n ( U ^"Wd)  = 0. 

m “ 1 

Let  ti.i  <j<  tj  and  Z^’  G S_.  Setting  Q = S.  in  Eq. 
(20),  we  have  {Xi^'>, . . . , X'„"’}  C (X^i  + K).  Setting 
g = S+  we  have  further  {X'^'’,  . . . X]''-'^} 

E (X^T  + K)  . That  is,  X^'i  ^ 

Let  r,-_i  < j < tj  and  G So . Substituting  g = 5o  into 
Eq.  (20)  we  get  G cone  (So),  j + 1 ^ p ^ n, 

together  with  So  = -So  and  So  H S+  = 0 this  yields 
{X„^'>,  ...,X‘„">}n(X+X^'')  = 0.  For  g = S\S^  and 


p G (1 , . . . , j-1  , Ui } we  have  X*^'  ^ 

{K  + X‘^').  Similarly,  for  p G {fi,r2,  • • • , Ui } we  have 

{Xi'>>,XM  ...,Xi''-'>}CXW  + X. 

Thus  in  this  case  X^}  G 

In  the  same  manner,  Z^^*  G S+\S_  implies  X'^'j}  G 

Summarizing,  if  Eq.  (18)  holds  then 
^'■>(^d  = (X*;-' : U,  <j^  tj,  Z“  G S^\S.},  whence 
Eq.  (19). 

Now  from  Eq.  (13)  and  Lemma  5 (recall  that  t = tk) 

P{(Z'‘’, . . . , Z''‘’)  G sis,  efi''’  > fif*; 
i — , tk  }—>p”‘{sis). 

Recalling  the  definitions  of  si  and  sis  we  get 
lim  inf  P{ = v,- , 7^  =f,- ; i^k}> 

n~^oo 

¥{V^'^  = Vi,P'^^tr,  i^kyS. 

Take  5 = 5(vi, . . . , v^;  ti, . . . , r*)  with 

.S5(vi, . . . ,Vk-,t\, . . . ,tk)  = (i  and  choose  a diagonal 
subsequence  of  the  values  of  n to  get  the  convergence  of 
the  probabilities  in  the  left-hand  side.  Recalling  that 
probabilities  sum  to  one,  we  derive  Eq.  (15)  by  setting 
/3^0  □ 

Convergence  in  mean  does  not  require  additional 
restrictions,  as  shown  next. 

Theorem  3.  Under  the  assumptions  of  Theorem  2 

Proof.  It  is  sufficient  to  consider  only  the  case  k’=  \. 
Denote  by  and  the  indicator  functionsof  the 
events  {X>'>  E ,2^'>(^d}  and  {Z''^  E S\S_,  1^7’"}, 
respectively.  Clearly, 

Vd"  = /'d+  . . . /'d,  + . . . 

By  an  argument  similar  to  that  used  in  Theorem  2 we 
show  that 

(/‘d,  ^ m=l,2, ...  (21) 

Choose  y„  satisfying  np„^A  , where  A > 0,  = 

The  random  variable  A = #({Xi, . . . , X„ } 
n ByJ  has  binomial  distribution  with  parameters  {n,p„). 
By  Eq.  (6), 

EC/d""^'’  + . . . + /’d)  = E(^'\^„)  r\ByJ^  e-'V. 

(22) 
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Fix  m and  write  the  expectation  as 

z^'  e s,\5_,  #({z">, . . . , z^^"}  n (s\s_))  = /-i } + 

Ev?  = E(/'i'  + . . . /'„">)  = E(/'''  + . . . + 

p{{z'“, . . . , z^^"}  n s,  = 0,  z'^'  G Si, 

E(/r"+.../':')l(;v.„}+E(/r 

#({z‘'', . . . , z^^"}  n (s\s_))  = / } = 

E(/r"  + ... 

p {#({z"', ...,z'>-"}n  So)  =i-i, 

The  first  term  converges  by  Eq.  (21): 

#({z"', . . . , z^"}  n (s_\s,))  = j-i}x 

E(/'l'+  . . . /'"')^E(/'"  + . . . The  first  and  the 

third  terms  are  estimated  by  Eq.  (22)  as 

PfZ^i  G S,\S_}  + P{#({Z">, . . . , Z^'^}  n So)  = i, 

#({Z'",.. .,  Z*^'*}  n (S„\S,))  = j-i-l  }P{Z^^  G Si}  = 

E(/r"  + ...  /'„"') i(/v^„>+E(/r"  + ...  /';') i(iv<„} 

^ E(/f^''  + . . . +/'„"')  ^ e-'V. 

{]l\)po'{p--p.r(p.-p.)  + (p--p.r-'p.. 

Since  N is  binomially  distributed,  we  have  for  the 
second  term 

where  / > 0,y  > 1 and  j ^ i.  Summing  over  j we  arrive 
at  the  limiting  distribution  of  points  in  the  ^th  layer: 

E(/r"''  + .../'i'')  i<^.„><E(yv  1<^^„>)-A 

P{t/(*)  = 0>  = — , 

l-p_  + Pi 

AP{A^a  — w-1}, 

where  is  a Poisson  random  variable  with  parameter 

A.  Selecting  A and  then  m sufficiently  large,  we  prove 
lim  sup„_»  EV‘^’  ^ EV<*'. 

The  inverse  inequality  involving  lim  inf  follows  from 
the  convergence  in  distribution  □ 

It  is  not  hard  to  find  the  limiting  distributions  of  the 
V*‘’’s.  Note  first  that 

pM/U:)  Po  (J~P-)P+  • 1 ^ 

If  the  cone  satisfies  K H -K  = 0 (or,  more  generally 
p(S+  n S_)  = 0 ) then  Pi  = 0 and  is  geometrically 

distributed.  A little  additional  work  is  needed  to  find  the 
expectation: 

(V"*  T")  (V*^*  J’(2)_J’(1)^ 

EV'*'  = . (23) 

7+ 

are  all  iid,  therefore  it  is  sufficient  to  consider  only  the 
first  pair.  Clearly  7"’  is  geometrically  distributed  with 
parameter  p(S+).  The  probability  law  of  V"’  is  found 
from  the  following  scheme:  throw  down  the  iid  points 
Z''',  ...  in  5 according  to  the  probability  law  p 

until  the  first  point  falls  into  S+  , then  count  all  the  points 
falling  into  5\S-.  To  make  this  precise  denote 

Example.  Assume  that  the  radial  tail  is  slowly  vary- 
ing and  p is  the  uniform  spherical  measure. 

For  ^ = Rf  we  have  p+  = p_  = 2^,  and  Eq.  (23)  yield 
EV‘^’  ^ 2“^  -1 . In  particular,  the  mean  number  of  Pareto 
points  in  two  dimensions  converges  to  3. 

Taking  the  complement  to  the  negative  orthant,  we 
have  p+  = p = 1-2“^  and  EV‘^*  — > (2*^-1 ) '.  In  two  dimen- 
sions, the  probability  that  the  sample  has  the  double 

p„  = p(So),  p±  = p(5i),  p_  = p(5-),  p^  = p(S+) 

maximum  tends  to  1/3. 

Appearing  of  inverse  numbers  in  the  above  example  is 
a general  phenomenon.We  write  further  V‘i’(^)  to 

(thus  po  + P-  + p+-p±  = 1 ).  The  joint  distribution  of  V"’ 
and  7*”  is  this: 

emphasize  the  dependence  on  the  cone. 

Theorem  4.  Under  the  assumptions  of  Theorem  2 

P{V‘"  = /,7“>  =j}  = P{{Z"', . . . , n 5+  =0, 

lim  EV'i*(/:)EV‘i>(-A:")=  1, 

G 5^,#({Z"',  ...,Z^'}n  (5\S_))  = / } = 

provided  one  of  the  numbers  p+  or  l-p_  is  positive. 

Proof.  This  follows  from  Theorem  3,  Eq.  (23),  and  the 
formulas  pf-K")  = l-p+(K),  p+i-K")  = l-pfK)  □ 

P{{z"', . . .,z^"}  ns,  =0, 
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Remark.  Given  a binary  relation,  say  S?,  on  a sam- 
pling space,  and  a random  sample  there  are  two 
natural  ways  to  define  the  “A-th  extremes”  of  : (1) 
sample  elements  X,  which  are  in  fk  with  all  other  sample 
elements  with  the  exception  of  some  k-\  points;  or 
(2)  the  elements  X,  such  that  there  are  exactly  k-\ 
sample  points  which  are  in  the  relation  with  X,  . In  the 
theory  of  partially  ordered  sets  extremes  (k  = 1 ) of  the 
first  type  are  called  the  greatest  points,  of  the  second 
type — maximal  [4].  This  is  best  illustrated  by  the 
natural  partial  order  of  R*^:  total  maximum  is  the  great- 
est point,  while  Pareto  set  consists  of  maximal  points.  If 
the  binary  relation  is  generated  by  a cone  X,  as 
mentioned  in  Introduction,  then  the  X -extremes  are 
maximal  points,  while  the  -X*^  -extremes  are  the  greatest 
points  w.r.t.  &t.  Baryshnikov  [3]  has  proved  that  the 
asymptotic  upper  bound  for  the  product  of  expectations 
of  the  numbers  of  the  extremes  of  both  types  is  at  most 
1,  for  any  fixed  S?  and  k.  Theorem  4 shows  that  this 
bound  is  sharp. 

Remark.  Normal  multivariate  distributions  can  be 
viewed  as  the  case  of  fast  decreasing  radial  tails,  a = oo. 
The  mean  number  of  conical  extremes  demonstrates 
typically  the  following  behavior:  for  any  k,  EV‘ 
infinitely  grows  if  X is  contained  in  a half-space,  and 
tends  to  zero  if  X contains  a half-space  [10,1 1]. 
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